THE AVERAGED CONTROL SYSTEM 
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Abstract. For control systems that either have an explicit periodic dependence on time or have 
periodic solutions and small controls, we define an average control system that takes into account 
all possible variations of the control, and prove that its solutions approximate all solutions of the 
oscillating system as oscillations go faster. 

The dimension of its velocity set is characterized geometrically. When it is maximum the average 
system defines a Finsler metric, unfortunately not twice differentiable in general. For minimum time 
control, this average system allows one to give a rigorous proof that averaging the Hamiltonian given 
by the maximum principle is a valid approximation. 
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1. Introduction. We consider eitlier a "fast-oscillating control system" ([T]): 

X = UiXi{-,x) -\ \-UmXrn{-,x) , ||w||<l, 

e £ 
where all Xis are 27r-periodic with respect to i/e, or a "Kepler control system" (|T7| : 

i = hiO + VlhiO + ■■■+ VmfraiO , \\v\\ < 8 

where all solutions of ^ = /o(^) are periodic. 

Averaging techniques for conservative — periodic or not — ordinary differential 
equations (ODEs) date back at least to H. Poincare; see |21 §52] or [^3] for recent 
expositions. Roughly speaking, on a fixed interval, the solutions of i = F{t/e,x) 
differ from these of i = Fi^) by a term of order e, with F the average of F with 
respect to its first argument. 

If u or v above is assigned to be a fixed function of state and time (or computed 
from additional state variables as in w = a{p,x), p = g{p,x)), then these techniques 
for ODEs can be applied to give an approximation at first order with respect to small e 
of the movement of the slow variables. Averaging is usually used in this way in control 
theory: in vibrational control [12], fast oscillating controls are designed and averaging 
techniques allows analysis and proof of stability; in the same way, it solves stability 
and path planning questions in control of mechanical systems, see for instance [8]; 
in |12[ §5], high frequency control is used to approach a non-flat system by a flat 
one; one may also mention many applications to control |21[ [T51 [2U] of the work [T7] 
that mimics Lie brackets by highly oscillatory controls along the original vector fields. 
A common feature to these references is that the use of oscillations "creates" new 
independent controls used for design. The use of averaging in optimal control of 
oscillating systems [101 [T^ [HI [7] is similar in spirit to the above, but closer to the 
framework of this paper because oscillations are present in the system instead of being 
introduced by the control. Very interesting results are obtained applying averaging to 
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the Hamiltoniaii equations arising from Pontryagin Maximum principle. For instance, 
in [7], the authors have studied in this way the problem of minimal energy transfer 
between two elliptic orbits; extremals are the same as those giving the geodesies 
of a Riemannian metric. Again, averaging introduces "new independent controls": 
Riemannian geodesies are minimizers of a problem where all velocity directions are 
allowed whereas the velocity set of the original system at each point had positive 
codimension. The same averaging computation may be applied to the Hamiltonian 
differential equation obtained for minimum time, but, since this differential equation 
is discontinuous, there is no theoretical justification for averaging in that case. 

Our contribution is to introduce a different way of averaging that takes into 
account all possible variations of the control — hence the control strategy can be 
decided after performing averaging — and to prove that it has satisfying regularity 
properties and is a good first order approximation of the above systems as e — > 0. 
This gives, as a side result, a justification of the use of averaging for minimum time 
in [T3l I14|. This procedure also "creates new independent control", i.e. increases the 
dimension of the velocity set, that we characterize in terms of the original vector fields. 
When this dimension is maximum, the average system defines a Finsler metric [3] on 
the manifold, whose geodesies are the limits of minimum time trajectories for the 
original systems as e — > 0. This Finsler metric is in general not twice differentiable 
(hence it is not a Finsler metric in the sense of [3], indeed); we however prove that, 
at least in the less degenerate case, the Hamiltonian system governing extremals, 
although it is not locally Lipschitz, generates a flow on the cotangent bundle. Low 
thrust planar orbit transfer belongs to this less degenerate case. 

The average control system may be used for other purposes than optimal control, 
for instance [4 designs a Lyapunov function for feedback control in the average system 
and uses it for the oscillating systems; indeed the present work was developed out 
of comparing feedback control based on a priori chosen Lyapunov functions with 
minimum time control for low thrust orbital transfer. 

Preliminary versions of this work can be found in [5J |3]. 

Organization of the paper. The construction and results are developed for "fast- 
oscillating control system" in ^and then transferred in ^to "Kepler control systems", 
and applied to minimum time orbit transfer in the planar 2-body problem in ^ One 
proof is postponed until the appendix, that comes after a short conclusion (^. 

2. Notations and conventions. 

2.1. M is a smooth connected manifold of dimension n; its tangent and cotangent 
bundles are denoted by TM and T*M. One may assume for simplicity M = M", 
TM = R" X M", T*M = R" X (M")*, and, for xG M, T^M = R", T*M = (R")*. 

For V G TxM, p G T*M (or any v,p taken in a vector space and its dual), we 
denote by (p, v) (rather than p{v)) their duality product. 

2.2. If £^ is a subset of a vector space V, then E-^ is its annihilator, the vector 
subspace of its dual V* made of all p's such that {p, v) — for all v in E. 

2.3. We assume that AI is endowed with an arbitrary Riemannian distance d. If 
M = R", just choose the canonical Euclidean distance. 

In local coordinates, ||.|| and (. |.) stand for the canonical Euclidean norm and 
scalar product. On a compact coordinate chart, ki\\x — y\\ < d{x,y) < k2\\x — y\\ for 
some positive fci, k2 (Lipschitz equivalence). We also denote operator norms by ||.||. 

2.4. is R/27rZ. For 9 in (an angle), we denote by fi{0) the unique real 
number in [0, 2tt) such that ij.(9) = 9 mod 2tt. For a real number s e R, we denote 
the angle it represents by s mod27r; it belongs to the quotient . 
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Maps S*^ i? (arbitrary set) are identified with 27r-periodic maps R ^ E. For 
instance, if / is such a map E and t € R, we write /(r) instead of /(r mod 27r); 

the average of / is denoted by ^ J^J{e)d9, or ^ J^^f{0)de, or ^ /^^^i /(6')d6'; one 
identifies LP(S'\R™) with the subset of LP(M,R™) made of 27r-periodic functions. 

2.5. The Euchdean norm in M™ or (M™)* is denoted by ||.||, and the bah of 
radius one centered at the origin by B™. We view an element of R"* as m x 1 
matrix (column) of real numbers and an element of (R*")* as a 1 x m matrix (line); 
transposition, denoted sends M™ to (M™)* and vice- versa. 

3. Fast oscillating control systems. We call fast oscillating control system 
on M a family of non-autonomous systems, linear in the control u G R™: 

t ™ t 

X ^ g{-,x)u ^y^gi{-,x)ui ,\\u\\ < 1 (1) 

indexed by a positive number e. Each Qi is a smooth "periodic time- varying" vector 
field: g, € C°°(S'i xM, TM). An admissible control is a measurable m(.) : [0,r] 
for some T > 0. For a given control u{.) and initial condition x(0), there is a unique 
solution x{.), defined either on [0,T] or only on a maximal interval [0,T'), T' < T. 

Remark 3.1. Apart from being a notation defined by the double equality in ([1]), 
Q{0, x) defines a linear map M™ — > T^M that sends (ui, . . . , Um)^ to X^I^li Qi{'^Tx)ui. 

3.1. Average control system of fast oscillating control systems. Define 

the map ^ : A/ x L°°([0, 27r],R") ^ TM by 



27r 



g[x,U) - / g(0,x)'U(f?)d0 . (2) 



It allows one to define, for all x G M, the subset £,{x) C TxM by: 

E{x) = {^(a;,!!), lie L°°([0,2^],R"), ||U||oo < l} C T^M , (3) 

and the average control system of ([T]) as follows; we explain in the following section 
to what extent it is the limit of ([1]) as e — > 0. 

Definition 3.2. The average control system of ^ is the differential inclusion 

X e £(x). (4) 

A solution of ^ is an absolutely continuous x{.) : [0,T] — > M such thatx(t) G E{x(t)) 
for almost all t. 

Proposition 3.3. For all x in M , £(a;) is convex, compact and symmetric with 
respect to the origin. 

Proof. It is closed, convex and symmetric because it is the image of the unit ball 
of L°° (S'^,R™) by a linear map; it is compact because g{x, .) is bounded on S*^. □ 

Further characterizations of £(a;) use the map H : T*M — [0, -l-oo) defined by 



2tt 



H{x,p) = ^l \\{p,g{9,x)}\\de (5) 



where the Euclidean norm is used according to §2.51 and, for each {6,x), 

{p,g{e,x)) = {{p,g,{e,x)),...,{p,g^{e,x))) e (r™)*. (6) 
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Proposition 3.4. For all {x,p) e T*M, one has, with H defined in ©, 

£(x) = |i;eT,M, sup {p,v) < l|, (7) 

H{x,p)<l 

H{x,p) = sup {p,v) = sup {p,g{x,U)) = (p,^(x,lt*.^)) , (8) 

veE{x) ueL°°(si ,»"•), ||ix|U<i 

. 1 s fo if {p,g(e,x)) =0, 

I ||(p,e(e,2:))|| y \Pi t'l'^j 2;;; 7^ u . 

Proof. The last equality in ([5]) is a straightforward maximization, the second 
one comes from the definition ([3]) of £(x) and a simple computation yields H{x,p) = 
(p, U*^)); this proves dH). Being closed and convex, E{x) is the intersection of 
all its supporting half-spaces jMl Corollary 1.3.5]; according to ([8]), this yields the 
following relation, equivalent to ([7]): E{x) = npiET*j\/ i"^ ^ T^jAf , {p,v) < H{x,p)}. □ 

A convenient characterization of solutions of (|4]) . According to Definition 13. 2i a 
solution x{.) is such that, for almost all t, there is Vi{t) e i°°([0, 27r], K"*) such that 
x{t) = g{x{t),Vi{t)); the map {t,d) 'U.{t){6) is measurable with respect to 6 but 
bears no regularity with respect to t. It turns out that it may always be chosen jointly 
measurable with respect to {t, 9) according to the following "measurable selection" 
result: 

Proposition 3.5. A map x : [0,T] -> M" is a solution of the differential in- 
clusion (HJ) if and only if there exists u G L°°([0,r] x 5'^,R™), |lw|ioo < 1 such that 



1 f^^ 

x{t) = — / g{e,x{t))u{t,9)de (lo) 

for almost all t in [0,r]. 

Proof. After possibly partitioning [0,T] into intervals where x{t) remains in the 
same coordinate chart, we work in coordinates and use a Euclidean norm when useful. 

Sufficiency is clear: from Fubini theorem, 9 i— )■ u(t, 9) is measurable for almost 
all t, hence x{.) is a solution of ([4]). Conversely, let a:(.) be a solution of (|4]): x{.) is 
measurable and, for almost all i, there exists ut G L°°{S^,W"'), \\ut\\oo < 1 such that 

x{t)^g{x{t),ut) = — j g[si,x(t))ut{si)dsi. (11) 

Let (j) : L°° ([0,T] x S-^M'") ^ ([0,T],M") be the linear map defined by 

1 /-^^ 

C^{u){t) = g{x{t),u{t, .)) = ^ Gisi,x{t)) U{t, Sl) dSl 

and J the image by of the unit ball of L°°([0,T] x 5\R"). Since, by ([TT]), x{.) is 
essentially bounded, it is in ([0, T] , R"); since J is closed and convex in that Hilbert 
space, the distance from i; to 3 is reached for a unique element ^ £ 3: 

^^cj){u), «GL°°([0,T]x5\M") , \\u\\l^<1. 

Let us prove by contradiction that ^ ~ x, i.e. x{.) G J; this will end the proof. 
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lix^^, one has, for all u in the unit baU of L°° ([0, T] x S^,R"^) , 

{x-^\^iu)-^iu))^, < (12) 

with equality only if <p{u) — (j){u). Define u by u{t^s) ~ ^/./ -/ s>jT (s) with 

U;^ defined by ®; clearly, u is in the unit baU of L°° ([0,r] x S'i,]R'"), and, for all 
(t,s) € [0,r] X 51 and ah u e R", 

||u||<l ^ (i(t)-C(t))^e(s,a;(t)) (u(i,s)-u)>0, (13) 

hence {x{t) — ^{t))^ G{si,x{t)) {u{t, si) — u{t, si)) is non-negative for almost all (i, si) 
and, since it is the integrand of the left-hand side of (|12p . it must be zero; hence 
^ — (t){u) — (j){u) and ^{t) = G{x{t),u{t, .)) for almost all t. 
In (jlip. ut satisfies ||Mt(si)|| < 1 for almost all si, hence, according to ([T3l) . 

(i(t)-C»)^e(si,x(t)) (^2(t,si)-wt(si)) >0. 

Since i;(i) = 5(x(t),'St), ^(i) = tj(a;(f), M(t, .)), the integration with respect to the 
variable si yields — ||i;(i) — ^ for almost all t; this contradicts x £,■ □ 

Remark 3. 6. The differential inclusion @ is equivalent to the "control system" 

x^g{x,U), UgL°°(51,R"), ||U||oo < 1 

where, by Proposition l3.5[ admissible controls are maps t h- > \i{t) such that u: {t, 9) 1-^ 
'U.{t){9) is measurable with respect to {t, 9). Since this "control" is infinite dimensional, 
and we could not find a representation of the type x = f{x,v), v € U G W, r finite, 
we stay with the differential inclusion with £,{x) described by ([SJ and 

3.2. Convergence theorem. The following result relates solutions of the fast 
oscillating systems as e tends to zero to solutions of the average system. To our 
knowledge, this kind of theorem where the control is not chosen prior to averaging 
has never been stated in the literature. 

Theorem 3.7 (Convergence for fast-oscillating control systems). 

1. Let Xo{.) : [0,T] — >■ M be an arbitrary solution of Q. There exist a family 
of measurable functions Ue{-) '■ [0,T] — > B"^ , indexed by s > 0, and positive constants 
c, £0, such that, calling Xe{.) the solution of ([T]) with control u = Ue{t) and initial 
condition Xe(0) = a;o(0), one has: Xg{.) is defined on [0,T] for all e smaller than 
£0 and converges to xo{-) as £ — > 0, with an error of uniform order e: 

d{x^{t),XQ{t)) < ce, [0,r], 0<£<£o. (14) 

2. Let K 6e a compact subset of M , (£„)„gN 0, decreasing sequence of positive 
real numbers converging to zero, and, for each n, a;„(.) : [0,r] — > K a solution of (J^l 
with e ~ En and for some control u — Un{t), Un{-) G L°°([0, T], R™), ||u„(.)|joo < 1- 
Then the sequence (3^n(-))„gN '^^ compact for the topology of uniform convergence on 
[0,T] and any accumulation point is a solution of the average system Q. 

The statement is more complex than the one for ODEs, e.g. [2,, §52. C], due to 
underdetermination (choice of control in ([T]), multi- valued right-hand side in (|3])). 

Informally, 'tJ]' states that any solution of the average system is the limit of 
solutions of fast oscillating systems with well chosen controls and 'l2]' states that, 
conversely, any limit of solutions of oscillating systems, with arbitrary controls, is a 
solution of the average control system. There is an estimate on the error in 'H]' but 
not in '12]' because some sequences may converge slower than others. 
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Remark 3.8. We could have considered systems that are affine instead of linear in 
the control by adding a drift vector field Qo{t/£, x) to ([T|). One then obtains an average 
control system where £ (x) is replaced by QQ{x) + t{x), withfJo(a^) — 9^ /q ^^0(^*7 2;)d0. 
Convergence does hold. We restrain to systems ([T]) because they naturally occur when 
dealing with a time- invariant system whose drift has periodic solutions, see §4] below. 



In the proof of Theorem 13.71 the following technical lemma is needed. 
Lemma 3.9. Let e > and a < b be real numbers and u : [a — 27r£, 6] x 
be measurable, one has the following identity (see §^2.4\ for the notation IJ,{.)) : 



g{0,x{s))u{s,e)d0ds^ g{-,x{s))u{s + en{e),-)d0ds + A, (15) 

<s<h a<s<b 



With A,^ II g{-,x{s + £n{0)))u{s + efiie),-)d9ds 

jpa 6 £ 



gi-,xis + en{e))) u{s + e^i{e), -)deds 



a<s<6 



g{-,x{s + £^i{e)))-g(-,x{s)) u{s + £^l{e),-)deds (i6) 



£ 



£ 



and the setT^ defined by T^ — {{s,d), OeS^, a~£^{6) < s < a} and accordingly. 

Proof. Thanks to the change of variables 6 = t/£ mod27r, s = r + e fJ,{(j>), with 
/i(6') as defined in ti2.4[ the left-hand side of ([T5|) is equal to 

g{-,x{T + £fi{(f))))u{T + £fi{(l)),-)dTd(j) . 

0t:O £ £ 

a— e(p<T<b— e i-L{(p) 

Keeping the names {s,6) instead of (r, 0), one gets ([T5|) . the correcting term 
coming from the modified domain of integration and argument of x. □ 
Proof of Theorem\3^ pomtUl Consider a solution xn : [0,T] -> Af" of g]). By 
Proposition [3H there exists uq e L°°([0,T] x S'i,K™), ||uo||oo < 1 satisfying For 
£>0, define Ue(.) e L°°([0, T], M") by (see 323] for notations): 

Mt)^^l uo{t + £KO),-)d0 , (17) 
27r Jo £ 

where uq is prolonged by zero outside [0, T]: uo{t + e /i(6'), |) = if t + e fi{d) > T. 
Let us prove that this construction of satisfies the two announced properties. 

Step 1. Let us first assume that M is an open subset o/M" and g is zero outside 
a compact subset of M . Then g{d^x) is a n x m matrix for all {0,x) and, denoting 
by ||.|| the Euclidean norm for vectors and the operator norm for matrices, there are 
global constants Lipt/ and sup^J such that, for all x, x' , 6 in M x M x S^, 

\\gi9,x)-gie,x')\\<{upg)\\x~x'\\, \\gid,x)\\ < (supg) . (is) 

Let 6 be a non-negative constant and consider, for each e > 0, a solution Xe{.) of ([T|) 
with control u — u^{t) and initial condition Xe(0) such that 

||x,(0)-xo(0)|| <6e. (19) 
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In fact 6 = in the theorem but we need a nonzero b in step 2. By definition, 
expanding Ite(s) as in (jl7p and using Lemma 13.91 one has 

x,{t)^xM + :^ f I e(-,a;e(s))uo(s-eM(^),-)d^ds, 
27r Jo Jo e e 

= a:e(0) + ^ I I g{e,Xe{s))uo{s,0)Aeds ^ (20) 



with Ae given by (fT6|) . that satisfies ||Ae|l < 47r^(Lip5) (l + TsupC/)e because, in 
particular, ||mo|| < 1, < ^tte and 

II {g{-^,xM) - Qi^^^eis + e ^i{e))))Ms + e ^i{e),'--)\\ <2n {Uj>g) (snj>g)e. 

Using (fT9|) . (|20|) the above bound on ||Ae|| and the relation 

xo(t) = xo(o) + -^ / / g{e.xo{s))uo{s-e)<ieAs, 

^TT Jo Jo 

yields 

||x,(i) -a;o(i)|| < (& + 27r(Lipg) (l + T sup g)) £ + (Lip 5) / ||x,(s) - a;o(s)|| ds 

Jo 

for all t in [0,r], and finally, using the classical Gronwall lemma, 

\\xS)-xo{t)\\ < [6 + 2^(Lipg) (I + Tsupg)] e^L'P^e (21) 

for all t in [0,T] and e in [0,£o]- This proves the theorem if M is an open subset of 
R" and g is zero outside a compact subset, with an explicit constant c corresponding 
to the distance d defined from the Euclidean norm and with Eq = +oo. 

Step 2. General case. Let x^{.) be the solution of ([ij with control u = Ug{t) 
defined in ()17p from uq and with initial condition Xe{Q) ~ a;o(0); it is not necessarily 
defined on [0, T] but may have a maximum interval of definition [0, T^) with < T. 
Let T e [0, T] be the supremum of the set of numbers re [0, T] such that, for some Eq 
and some c, that may depend on r, the solution Xe(.) is defined on [0,t] and satisfies 
d(a;e(i), xo(0)) < ce for all t G [0,t] and e G [0,eo]- Let us prove by contradiction 
that T = T. This will end the proof of Theorem EUl point □ 

Assume T <T, and let _ 

- O be a coordinate neighborhood of xo{T), 

- a > be such that 0<f-a<f + a<Tand xo{[f - a, T + a]) C O, 

- c > 0, £0 > be such that d{x^{t), xo{t)) < ce for all t e [0, T — a] and e £ [0, £o]- 
Taking £o possibly smaller, one also has Xe{T~a) £ O for £ < £o. Let K be a compact 
neighborhood of xo{[T — a,T + a]) contained in O, K' a compact neighborhood of K 
contained in O, and p : M — >■ [0, 1] a smooth map, zero outside K' and constant equal 
to 1 in K. Defining by gp{9,x) = p{x)g{9,x), let us apply Step 1 in coordinates 
in O, with C/^ instead of g and [T — a,T + a] instead of [0,T]. Call Xq (resp. x§, 
£ > 0) the solution of (jlOp (resp. of ([T]) with control u = Ue{t)), replacing g by C/p, 
with initial condition x^{T ~ a) — Xe{T ~ a), e > 0. One clearly has, as in 
||a;^(r — a) — a;o(r — q;)|| < be with b deduced from c via Lipschitz equivalence of the 
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distance d and the one induced by the norm in coordinates (see iJ2.3p : then Step 1 
provides > and, by (l?T|) . the inequaUty 

|K(t)-xg(t)|| < [b + 2Tr {Up g + 2a sup Gp)] e^^LipSpg^ (22) 

vahd for t E [T ~ a,T + a] and e E [0, Eq] . Possibly choosing a smaher £g, this impUes 
that Xe{[T — a,T + a]) C K for e < Eq; since Q that coincides with Gp in K, the 
conclusions holds for and Q as well as x^ and t/p if e is no larger than e'q. We 
have shown that, for all e < e'q, the solution Xg is defined on [0,r + a] and, with 
c' deduced from c on time time-interval [0, T — a] and from the explicit constant in 
(j22p and Lipschitz equivalence of d and the Euclidean distance on time time-interval 
[T — a,T + a], it satisfies d{x^{t) — XQ{t)) < c'e for t in [0, T + a]. This contradicts 
the definition of T. □ 

Proof of Theorem \3. 7\ point \M Since g is bounded on 5*^ x IK (one may cover 
K with a finite number of coordinate charts and define this bound in coordinates), 
the maps Xn(-) have a common Lipschitz constant and the sequence (a;„(.)) is equi- 
continuous, hence compact by Ascoli-Arzela Theorem: one may extract a uniformly 
convergent sub-sequence. Still denoting by (a;„(.))„gN such a converging sub-sequence 
and by x*{.) its (uniform) limit, we need to prove that this limit is a solution of Q. 

Define, for each n, m„ : [0, T] x S'^ ^ M™ by 

Un{t,e)=Un{Mt,ff)), (23) 

where Un{ ) G L°°([0, T], E™) is associated to Xn{ ) according to the assumption of 
the theorem and where the map /3n ■ [0, T] x S"^ — ?> M is defined by 

t-2TTEn< Pn{t,e)<t , ^"(^'^) =6> modulo 27r . (24) 

En 

Clearly u„ is in _L°°([0,r] x S'-'^,M'") and ||u„||oo < 1- Hence, after possibly extracting 
a sub-sequence, (u„) converges in the weak-* topology to some u* . Let us prove that, 
for almost all tE [0,T], 

1 f^^ 

ii) = Tr S{0, X* (t)) u* (t, 9)de . (25) 
27r Jo 

Let T E [0, T] be the supremum of the set of numbers r E [0, T] such that this is true 
for almost all t in [0, r], and let us prove by contradiction that T = T. 

Assume T < T, and let C be a coordinate neighborhood of xq{T) and a be such 
that 0<T ~ a<T + a<T and xq{[T ~ a,T + a]) cO. Uniform convergence implies 
Xni[T—a, T+a\) C O for n large enough and then, in coordinates, for t E [T—a, T+a], 

~ s 
Xn{t) - Xn{T ~ a) ^ g{ — ,a;„(s))u„(s)ds. (26) 

JT-a £n 

From (P^ . one has /J„(s -I- End, j-) = s, hence, from u„(s -l- e„6', ^) = u„(s) for 
all 6* e 5^, s £ M; using this in Lemma [3.91 one has 



27r 



/ g{e,Xn{s))unis,9)d0ds^ ^ 1 1 ^ ( — , (s) ) u„ (s) d(? ds + A,„ 

lT~a<s<t J JT-a<s<t 

ees^ ees^ 
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Since the integral in the right-hand side — whose integrand does not depend on 9 — 
is equal to the right-hand side of (pS)) . one gets, using uniform convergence of x„ to 
X* , weak convergence of u„ to u* and convergence of A^^^ to zero, 

x*{t)-x*{f-a) = ^[ [ g{9,x*{s))u*{s,e)d0ds , 

for t in [f - a,f + a], and finally that hold for almost all t in [0,f + a], thus 
contradicting the definition of T. □ 

3.3. Dimension of the velocity set £,{x). Recall that, for a convex subset C 
of a linear space, containing the origin, its linear hull is the smallest linear subspace 
that contains C, the interior of C in its linear hull is always nonempty, dim C is the 
dimension of this linear hull. 

Viewing ^~-(9,x) as a linear map R™ TxM (see Remark l3.ip . and S denoting 
a sum of linear subspaces of T^M , define the integer r{9, x) by: 

r(0,a;) =dim(^ Range 1^(0, a;)) . (27) 

It is also the rank of the collection of vectors ^^f{9, x) £ T^jM, 1 < i < to, j > 0. 

In the following proposition, and it is the sole place where this property is used, 
"system ((Ij is real analytic with respect to 6'" means that the vector fields Qi are real 
analytic with respect to 9 for fixed x (while being smooth with respect to (0, a;)). 

Proposition 3.10. 

1. The linear hull of £(a;) satisfies the following two properties for all x in M , 
where the inclusion Ii29\} is an equality if (Q]) is real analytic with respect to 9: 

Linear hull £ (a;) = ^ Range5(6', a;) , (28) 
Linear hull £ (a;) D ^ Range -^(6*, x) for all 9 e . (29) 

2. Ifr{9,x) — n for at least one 9 in , then £(a;) has a nonempty interior in 
TxM , i.e. dim£(a;) = n. 

3. If the system {!]) is real analytic with respect to 9, then r(9, x) does not depend 
on 9 and r{9,x) = dim£(a;). 

Proof. If p is in Ka.nge Q {9 , x) for all 9, then any v = Q{x,Vl) in £(a;) sat- 
isfies {p,v) = because {p,Q{9,x)Vi{9)) is identically zero on [0,27r]. Conversely, 
let p be in £(a;)^, and consider v — Q{x,Vip^) G £(a;); then {p,v) = implies 

{p,g{9,x)) = 0, i.e. p e Rangee;(6',x)^ for all 9. This proves the identity (j28p 
between subspaces of T^M through the following one on their annihilators in T*M: 
£(a:)-'- ~ Pleesi (Rs.nge G {9 , x)) . Similarly, and using the inclusion implies 

n^esi (Range e?((/),a;))^ C f]j^f^{Ra,nge ^{9,x))'^ . lip is in Range^((/), a;) for all 
0, differentiating {p , G{(j), x)) — with respect to cf) an arbitrary number of times yields 
{p ,d^Q /d(t)^{x, (f))) — 0, j & N; this proves, with = 6', the above inclusion, hence p9|. 

To prove the converse in the real analytic case, fix 6* e 5*^ and p e HjeN ^) ' 

and consider the real analytic mapping (M'")*, {p ,Q{(f),x)); the assump- 

tion on p implies that this map vanishes for (f> — 9, as well as its derivatives at all 
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orders, hence is identically zero: p e H^esi (^^^E^S{4>, x))'^ . This ends the proof of 
Point [TJ Point [5] is an easy consequence and Point [3] is classical. 

3.4. Further properties in the full rank case. We now assume that the 
mapping in ([T]) is such that the rank r{9, x) defined by (j27p is maximal: 

r{e, x) =n for all x in M and 9 in 5*^ (30) 

3.4.1. Controllability. This condition is strongly related to controllability of 
the linear approximation of ^ around equilibria, i.e. around the solutions where x 
constant and u identically zero. Indeed, the linear approximation of the time-varying 
nonlinear system (take e = 1 in ([T])): 

x = Q{t, x)u (31) 

around the equilibrium x = xi \& the time-varying linear system ^ = Q{t,x\)u\ ac- 
cording to [16l p. 614], it is "controllable with impulsive controls at any time" if and 
only if r(t, xi) = n for all t. If this is true at all points xi then all end-point mappings 
are submersions around zero controls; we shall need the following more precise result: 
Proposition 3.11. Assume that iT^) holds. 

1. For all xi S M and T>0, there exist a coordinate neighborhood W of xi 
(the ball B below refers to the Euclidean norm in these coordinates), positive constants 
q;o,C3, and, for all y G W, a smooth map Xy ■ 'S(j/,ao) — ^ L°°([0, T], K™) with 
Lipschitz constant c^, which is a right inverse of the end-point mapping of iSl]) on 
[0,T] starting from y, i.e. for all yf G B{y,ao), the control XyiUf) ■ [OjT'] — > M™ is 
such that the solution of x = Q{t, x)xy{y f){t) i a^(0) = y satisfies x(T) — yf. 

2. For all e > 0, the system (QP is fully controllable, i.e. there exists, for any 
e > and any two point Xa,Xi in M, a time T and a measurable control u : [0,T] — > 
_B™ such that the solution of |IJ) with x{0) — Xq satisfies x{T) = xi. 

Proof Let Ey : L°°([0, T], M™) M be the end-point mapping with starting 
point y. Condition (I30p implies that the derivative of E^-^ at the zero control has rank 
n; hence there exists an n-dimensional subspace V of L°°([0, T], R™) such that the 
restriction of E^-^, and hence of Ey for y close enough, to is a local diffeomorphism 
at zero; the Xy's are the local inverses of these local diffeomorphisms; they depend 
smoothly on y, hence the common ao and C3 in Point 1. 

This implies that the reachable set from any point at any positive or negative 
time contains a neighborhood of this point; a classical argument then tells us that 
the reachable set from a point xq is M , assumed to be connected, for it is both open 
(obvious) and closed (if x is in the closure of the reachable set, some points in the 
reachable set can be reached in negative time, hence x can be reached from xq). □ 

Let us now turn to the average system From H : T*M [0, +00) defined 
by ©, we define N : TM -J> [0, +00] by 

N(x,v) = max (p,v) . (32) 



Proposition 3.12. Assume the rank condition (3^. 

1. For all X G M , H(x, .) defines a norm on the cotangent space T*M, its dual 
norm on the tangent space T^M is N{x, .), and £,{x) is the unit ball for N{x, .), i.e. 
£(x) = {ve T^M, N{x,v) < 1}. 

2. System ^ is fully controllable, i.e. there exists, for any points xo,xi in M, 
a time T and a solution x{.) : [0,T] — M of ^ such that x{0) = xq, x{T) = xi. 
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Proof. From H{x,p) — implies {p,Q{d,x)) = for all 6 and, differentiating 
with respect to 6 and using pop . this implies p = 0; this makes p i— > H{x,p) a 
norm for other properties are straightforward. Hence N given by (|32p is finite for 
any {x^v) and it is, by definition, the dual norm of H{x,.); £,{x) is its mrit ball by 
^ in Proposition 13.41 To prove Point [2l take a continuously differentiable curve 
7 : [0,1] — ^ M such that 7(0) = xq and 7(1) = a;i and a : [0,T] — > [0,1] for some 
T > 0, differentiable, such that 

Jo ds 

(A'' and H are obviously continuous), then t a:(t) = 7(fT(t)) is a solution of such 
that a;(0) = xq and x{T) — xi. □ 

3.4.2. On the differentiability of if. It is clear that H, given by is as 
smooth as G on T*A/ \ Z (where the square root does not vanish) with 

Z - {{x,p) e T*M, 36 e S\ (p, Gie, x)) = 0} . (33) 

Unfortunately Z is not empty: it is, generically (away from {p — 0}) a 2n — m + 1 
dimensional submanifold of T*M. Let us give a first statement valid even on Z. We 
state it for : (x,p) H{x,p)'^, because H itself, homogeneous of degree 1 with 
respect to p, cannot be differentiable on {p — 0}, that coincides with {H{x,p) ~ 0} 
by Proposition 13.121 item [TJ 

Theorem 3.13. If condition i30\) holds, is continuously differentiable. 

We first prove a more generic property. In the following proposition, the notations 
are independent from the rest of the paper: 

Proposition 3.14. Let d be a positive integer, O'^ an open subset o/M^, V: 
5^ X O'^ — )• M™ a smooth map (C°°), Z the subset of O'^ where V vanishes for some 
9 and H: O*^— > [0, +00) the average of the norm ofV: 

H(x) = — ||v(e',x)||de' , z = {x eO''.39es^,\/{e,x)^o} . (34) 

27r Jo 

Assume: 

for all X in O'^ , the set {9 Cz , \/{9,X) = 0} has measure zero in . (35) 



Then H is continuously differentiable and, for all X , 

r-27r / y{9,X) 



1 r^"" f d\l 



Mo,x)\\ 



d9 . (36) 



Proof. The integral in ([36| is well defined (its integrand is bounded) and, by ((35|) 
and Lebesgue convergence theorem, it is continuous with respect to X and h. Let us 
prove that this dH is the derivative of H. Since V is smooth, one has 

\\y{9,X + h) - \/{9,X) - —i9,X).h\\ < k\\h\\\ (37) 

where ^{9,X) is smooth with respect to {9,X) and k is some local constant. Now, 
assuming \/{9,X) ^ 0, one has 

y{9,X) 



||V(0,X + /t)|| - ||V(0,X)|| = \ \l{9,X^h)-y{9,X) 

+ a{9,X, h) - 



yi9,x)\\ 
y{9,x + h) -y{9,x)\\^ 



\/i9,X)\\ + \\yi9,X + h)\\ 
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with \a{9,X,h)\ < 2. Hence, from ([57)) and (1551) . one has, for some local constant fc', 
\\H{X + h)-H{X)~dHiX).h\\ < fc^ A^^^^ ^ ||V(g,X + /^)-V(^,X)|| \ ^ 



- 27rJo V ||V(0,X)|| + ||V(0,X + /i)|| 

for \\h\\ small enough. For fixed X and ft, — > 0, the integrand in the right-hand side is 
bomided by 1 + |jft|| and converges to zero for 9 outside the set {6 £ S^, V(6', X) =0}: 
by ([35]) and Lebesgue convergence theorem, the right-hand side tends to zero. □ 

Proof of Theorem \3.1Sl We operate in coordinates for this is anyway a local 
property. Apply Proposition 13.141 with d — 2n, O'^ a neighborhood of a point where 
p ^ 0, X = {x,p) e M^" and \/{9,X) = {p,Q{0,x)): the rank condition implies that 
derivatives of all orders of the smooth map i-> \/{6,X) never vanish at the same 
point, so that its zeroes are isolated and the set {9 G , \/{0, X) = 0} is finite and 
a fortiori has measure zero; hence H is continuously difFerentiable outside {p — 0}. 
Since < H{x,p) < k\\p\\ for some local constant k the derivative of is zero at all 
points (x, 0) and, since p6p implies that the norm ofdH(x,p) at neighboring points 
where p ^ is bounded, the derivative of at these points tends to zero as p — 0. 

is therefore continuously difFerentiable everywhere. □ 

The map H fails in general to be twice difFerentiable on Z; indeed its first deriva- 
tive fails to be Lipschitz continuous, but we have the following estimate of its modulus 
of continuity, in the less degenerated case. 

Theorem 3.15. Assume that the rank condition Ii30\} holds and that 

(i) for {x,p) £ T* M , p 7^ 0, there is at most one 0(z such that (p, 0(0, x)) = 0, 
dQ 

and (p, -—(0,x)) does not vanish at the same point, 
o6 

(ii) for all {9,x) G S*^ x M, one has rankC/(0, x) = m, 

then any point {x,p) has a constant c and a coordinate neighborhood in T*M such 
that for all X and Y in R^", coordinates of points in the neighborhood, 

\\dH{Y) - dHiX)\\ <c\\X- Y\\ In ■ (38) 

In the left-hand side of p8l) . and PO)) below, ||.|| stands for the operator norm 
in coordinates, see ^2.31 As for Theorem I3.13[ we first prove a more generic result: 
Proposition 13.161 below, whose notations are these of Proposition 13.141 independent 
from the rest of the paper; its proof is in the Appendix to the paper. 

Proposition 3.16. Let V satisfy the following assumptions: 

(a) V(6l, X) ^0 for all {9, X) € S'^ x such that V(6', X) = 0, ^ 

(b) for any X E O'^ , there is at most one 6 such that \I{9,X) — Q^ > (39) 

(c) V and dM / 39 do not vanish simultaneously, ) 

and let X be in Z . There is a, neighborhood U of X in , and a constant K > 
such that, for all X, Y in U , 

dHiX)-dH{Y)\\<K\\X-Y\\ ln ||^ly|| ■ (40) 

Proof of Theorem \3.15[ Smoothness outside Z is obvious from the expression ([5]) 
of H; inequality (138p is a consequence of Proposition 13. 161 applied with d = 2n, X = 
{x,p) e ]R2«^ v(6', X) = (p, g{9, x)) and O'^ a neighborhood of a point of Z \ {p = 0}; 
it is clear that points (i) and (ii) imply the three conditions (|39|l . □ 
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Remark 3.17 (Finsler geometry). If was an least twice continuously differen- 
tiable, with a positive definite Hessian with respect to p, so would be iV^ (see (15^ ). 
and it would define a (reversible) Finsler metric [3J on AI . The lack of differentiability 
calls for further developments. 

3.5. Application to the minimum time problem. Fix two points xo,xi in 
M and consider the time optimal problem associated to (IlJ for e > 0: 

(v.), e>0: = ^(i/e, .(0).(t), u(i) E t e [0,T], 1 ^^^.^^ 

x(0) = xo, a;(T) = xi J 

and the time optimal problem associated to the average system: 

x{t) e £{x{t)), te[0,T], \ 

x{0) ^ xo, x{T) ^ X, / (^2) 

Call Ti;(x(j,xi) the minimum time for (Ve), £ > and To(a;o,a;i) the one for {Vq); 
when no confusion arises, we write Tg and Tq. 

Let us develop (|3T])-(|35]): concerning (|^T|) . is the infimum of the set of T's such 
that there is an admissible control u{.): [0,T]— and x{.): [0,T]— >-Af satisfying 
a;(0) = xo, x{T) = xi and x{t) = Q{t/e,x{t))u{t) for almost all t; Proposition 13.111 
point [5] implies that this set is nonempty, hence is finite. Concerning (jH)), Tq is 
the infimum of the set of T's such that there is x(.) : [0, T] M satisfying a;(0) — xq, 
x{T) = xi and x{t) G ?,{x{t)) for almost all t, Tq is finite from Proposition 13. 12[ point 
[2] A solution to {Ve) (resp. to [Vq)) is x{.),u{.) (resp. x(.)) as above with T — 
(resp. T = Tq). In general, the minimum T^ or Tq need not be reached, i.e. there 
need not be a solution. 

Lemma 3.18. Assume the rank condition I13(J\} . 

1. There is a neighborhood W of any xi and two constants ao > and C3 > 
such that, for all y in W, Ti;{y, xi) < 27re + C^dlxi, y). 

2. For any xq,Xi,Xi in M , one has T^{xq,Xi) <Ti;{xq,Xi) +Tg{xi,Xi) + 2716. 
Proof. Apply Proposition 13.111 point [1] with T — 2tt, using as a distance in W 

the Euclidean norm in some coordinates: for any two points y, y' in W such that 
\\y ^ y'W < c^o, there is a control defined on [0,27r], with L°° norm smaller than 
cslly — y'll that brings y' at time to y at time 2tt for system ((3T|) : rescaling time and 
control by e yields, if csWy — y'\\ < e, a control with L°° norm less than 1 that brings 
y' at time to ?/ at time 27re for system ([T]) and hence, by concatenating controls 
and using periodicity of for any positive integer k, a control with L°° norm less 
than 1 that brings y' at time to y at time 2fc7r£ for system ([T]) if c^Wy — y'\\ < fee. 
In other words, T^{y',y) < 27r(e + c^Wy — y'\\). Take y' — xi and 2ttc3/C3 the ratio 
between the Euclidean norm and the distance d; this proves point 1. Point 2 follows 
from using periodicity of Q and concatenating controls while inserting a zero control 
between time Ti;{y' , y) and the next multiple of 27r. 

Theorem 3.19 (limit of minimum time). Assume the rank condition VJO) . 

1. is hounded as e — !■ and limsupg_j.o T^ <Tq . 

2. //, for e > small enough, each (Ve) has a solution x^ : [0,Tg] M and 
there exists a compact K C M such that Xe{[0,Ts]) C K for all e > small enough, 
then all accumulation points of the compact family (xe(.))e>o in C^{[0,Tq], M) are 
solutions of (Vq) and — Tq . 

Proof. Consider a minimizing sequence for problem (Vq), i.e. solutions x'^ : 
[0, To + (3k] M of the average system ^ with (/?&) a sequence of positive numbers 
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that tends to zero and x*''(0) = Xq, x^{Tq + /J^.) = xi for all k. For each x*^, there is, 
according to Theorem 13. 7[ a family )£>o such that each x^{.) is a solution of 

^ with x^iQ) = xo and (i(a;^'(0> < 

ciE for all t in [0,To + /3fc]. In particular 
d(a;^(To + ^fc),x'=(t)i) < ci£. Now, from Lemma [XH] stated below, T^{x^;{Ta + 
f3k),xi) < (27r + ciC'3)e; hence, from the second point of that lemma (with x'^ = 
x^,{To + Pk) ), one has Te = Te(xo, xi) < To + + {An + ciCz)e and, letting k go to 
infinity, < To + (47r + ciC3)e; this implies Point [H 

In Point [5J first extend x^ on [0, T], with T an upperbound of T^, by taking 
Xgit) = xi for all t G [Te,r]; any sequence (xej.(.))fegN with limefc = is compact 
in C°([0, r], M); take a convergent subsequence such that T^^ also converges to some 
T*; from Theorem 13.71 the uniform limit is a solution of the average system and 
it goes through xq at time and xi at time T* , hence T* > Tq by definition of Tq. 
Since T* can be any accumulation point of (Tg) as e — )■ 0, this, together with Point [1] 
ends the proof of Point [21 □ 

Let us now write the Pontryagin Maximum Principle (PMP) [22] both for (Ve), 
e > and for (Vq) and see how they are related. 

The extremals of problem (Ve), e > 0, are absolutely continuous maps t i— >■ 
{x{t),p{t)) solution to 

p = i = -g^ with H,{t,p,x)^\\{p,g{t/e,x))\\, (43) 

whose right-hand side is discontinuous on = {{x,p,t), {p,Q(t/e,x)) — 0} (the 
"switching surface"), where it is in fact not defined. 

The extremals of (Vq) are absolutely continuous t h- !■ {x{t),p{t)) solution to 

with H given by (O. The right-hand sides are continuous according to Theorem 13. 131 

Theorem 3.20. // an absolutely continuous map t i— > x(t) defined on [0,T] is a 
solution of {Vs:), £ > Q, (resp. of (Vq)), then there exists t i— > p(t) defined on [0,T] 
such that 1 1— > {p{t),x{t)) is an extremal of (Ve), £ > (resp. of (Vq) ). 

Proof. Problem [Ve), £ > deals with a classical smooth control system; accord- 
ing to [251 E]) the pseudo-Hamiltonian is h{t^ x,p, u) = {p, G{t/e, x) u): an extremal is 
a curve on the co-tangent bundle solution, in local coordinates, of: 



p = -^{t,x,p,u*) = -{p,%u 
X = ^{t,x,p,u*) = Gu 



with u*{t) a control that maximizes the pseudo-Hamiltonian for almost all time; it is 
defined by u* = if (p, G{t/£, x)) 7^ 0; the maximized Hamiltonian He{t,p, x) — 

maxu h{t,x,p,u) is the one in (|43l) and (j45|) is then the differential equation (|43| . 
whose right-hand side is discontinuous at points where {p,G(t/e,x)) vanishes. 

Let us now turn to (Vq)- Since the set of admissible velocities is not a priori 
smooth with respect to the state variable we use a non-smooth version of the Pontrya- 
gin maximum principle for differential inclusions, that we recall for self-containedness: 
Theorem 9.1 in [llj Chapter 4]: if x € £-{x) is a locally Lipschitz differential 
inclusion and t 1— >■ x(t) is an absolutely continuous function defined on [0,T] 
solution to the problem (j42p . then there exists t 1— >■ p{t) defined on [0,T] such 
that {—p,x) G dcH{x,p) for almost allt G [0,T] with H{x,p) = max^g£(2,)(p, ij) 
anddcH the generalized gradient of H . 
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The set- valued map £(.) in ^ is indeed locally Lipschitz: in local coordinates , for 
xi, X2 in M", denoting by S the Hausdorff distance between two sets, one has: 



(5 (£(a;i), £(a;2)) = max < sup inf — W2II, sup inf ||wi 

= max min \\Qixi,Ui) - Q{x2,U2)\ 

||Ul||oo<l ||U2||oo<l 

According to ([5]), the Haniiltonian H defined in the above quoted theorem coincides 
with the map H defined in ([S])- □ 

Remark 3.21. This result and Theorem l3.19l have two interpretations: 

1. They prove that the operations of averaging and computing the Hamiltonian for 
the minimum time problem commute. Indeed, the Hamiltonian H was obtained by 
applying the maximum principle to problem (|42l) . i.e. minimum time for the average 
system ([4]), but it also the average of the one in (|43p with respect to the fast variable. 

2. They prove indirectly an averaging result for the minimum time control problem 
(|42p ; the averaging techniques in |10| do not apply to minimum time for they require 
smoothness of the Hamiltonian, while averaging is used in ITB] for minimum time 
with only partial theoretical justifications but numerical evidence of efficiency. 



Let us now focus on the differential equations (|44p that govern the extremals of 
(Po)- It is of great importance to know whether it defines a Hamiltonian flow on T*M, 
i.e. whether solutions trough all initial conditions are unique or not. Its right-hand 
side is continuous because, from Theorem 13. 131 H is continuously differentiable; this 
ensures existence of solutions. We saw that H is smooth {C°°) on T*M\Z (see (1551) ). 
hence solutions through points outside Z are always unique. The following result 
gives uniqueness of solutions even on Z in the less degenerated case possible. 

Theorem 3.22 (Hamiltonian flow for (T^o))- Assume that the rank condition 
^30fl holds, as well as conditions (i) and (ii) in Theorem \3.15i Then the differential 
equation (|44p has a unique solution from any initial condition. 

Proof. For an autonomous ODE i = f{z) in a finite dimensional space, where / 
satisfies ||/(zi) — /(z2)|| < t^(||^;i — ^2!!) with a; : [0, +03) — ^ [0, +00) non-decreasing, 
Kamke uniqueness Theorem [15, chap. HI, Th. 6.1] states that uniqueness of solutions 
holds if Jq ■;^li^ = +oo for arbitrarily small a > 0. From Theorem l3.15[ we are in this 

case with u>{u) = culn{l/u), and J — ~c Inln^. □ 

The sufficient condition given by this theorem is certainly not the best possible, 
but turns out to be applicable to the control of orbit transfer with low thrust, see |5l 
Point (ii) is very mild and only states that the control vector fields are linearly inde- 
pendent. Point (i) is more artificial: the fact that {p,Q{6,x)) = has at most one 
solution 9 has to be checked by hand, while the fact that dQ/dO does not vanish at 
the same time, which is true for the Kepler problem and used in [S] to show that 
the switchings in (|43)) are always "7r-singularities", i.e. the control u* switches to its 
opposite, is equivalent to the rank condition 

rank|g(0,a;),^(0,a;)| = dim(Rangeg(0, x) + Range ^(0, x)) = n . (46) 



Proving existence of a fiow for (j44p in more general situations is obviously an 
interesting program to be pursued. 
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4. Kepler control systems. We call Kepler control system with small control 
a family of control system on x M of the form 

where G and g can be viewed, with the same convention is in ([T]), as n x to and 1 x to 
matrices smoothly depending on {0, x) and w is a smooth function x M ^ M. that 
remains larger than a strictly positive constant: 

ij{9, x)>k^>0 V(6i, x) e X M . (48) 

In fact, this is an afhne control system on x M 

m 

^ = /o(0+E"'/'(^) (49) 

with ^ = (0, x), /o = w ^ and, for 1 < i < m. the smooth vector field fi is represented 
by the i'^ column of the matrix notations G and g. 

In if all solutions of ^ = /o(C) periodic but ^ is not a priori decomposed 
as a product, additional conditions are needed for the orbits to induce a nice foliation 
that splits the state manifold into a product M x S^. 

4.1. Relation with fast oscillating systems. Let us define a new time A in 
order to transform the first equation in (|47| into a "fast oscillating system" of the 
type ([T]). For a solution ti-^ {9{t),x{t)), let 0(t) be the cumulated angle i.e. Q{.) is 
continuous [0,T] -J> R with e{t) = 0{t) mod 2tt for all i and 9(0) S [0,27r), and let 

A = 3?(t) ^ £ (e(t) - 9(0)) . (50) 

Taking £o small enough that \uj{6,x) + e g{9,x)u\ > for a: £ K, ||u|| < 1, 

£ < £o, one has dJl/dt > ek^/2 hence 31 is strictly increasing and one-to-one, and 

k 

— £ t < JKt) < kuj et with fc^ sup co + Eq sup H^H . (51) 

2 S^xK S^xK 

Define, from'i;(.), the control A i— > u(A)= v{'3l^^{X))/e; then a;(A) = a;(3?^^(A)) satisfies 



dx _ G{0Q + j,x)u 

dX ~ uj{9o + ^,x) + eg{eo + ^,x)u 



i^Oo,e) — - V. , X ~^ , . A ^ < 1 • (52) 



Except for the term egu in the denominator, this is a fast oscillating system ([T]) with 
a suitable definition of Q. We will now apply the concepts and results of ^ 



4.2. Average control system. The definition uses uj defined by 

1 1 r^'^ d9 



uj{x) 2tt Jq uj{9, x) 



(53) 



Definition 4.1 (Average control system of Kepler control systems) . The average 
control .system of the Kepler control system is the differential inclusion 



X G E-{x) 



(54) 
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with £ defined by (0) using G : M x L°°(5\R") TM defined by 

Qix^U) = Zo{x)^r^^U{e)de (55) 



instead of ([2]). Solutions are defined as in Definition \3.Sl 

Remark 4-2. This Definition 13.21 applied to (15^ . which is equivalent to (H7| via 
time changes, except: 

(i) the term sgu in the denominator of ()52p has been discarded, 

(ii) the right-hand side has been multiplied by iJ{x) so that the "time" of the 
average system is equivalent to e times the "time" of the control system (j47|) as e ^ 0. 

4.3. Convergence Theorem. The counterpart of Theorem 13.71 is: 
Theorem 4.3 (Convergence for Kepler control systems). 

1. Let Xo{-) '■ [0,T] M be en arbitrary solution of (j54p and 9^ G . There 
exist a family of measurable functions u^{.) : [0,T] — > B™, indexed by s > 0, and 
positive constants c, eo; such that, if t (6s{t), Xe{t)) is the solution of (|47p with 
control u = Ue{t) and initial condition (0^(0), a;e(0)) = (0'^,xo(O)), it is defined on 
[0,r/e] for e smaller than Eq and 

T 

d{xe{t) , xo{et)) < ce , te[0,—], < e < eq , (56) 

thus T Xeir/e) converges uniformly on [0,T] to r h- )■ xo{t) when e tends to zero. 

2. Let K 6e a compact subset of M , (en)neN a decreasing sequence of positive 
real numbers converging to zero, and (0„(.), a;„(.)) : [0,r/e„] — > S*"*^ x K a solu- 
tion of system \4-T^ for each n, with e — Sn o,nd some control u = Un{t), «„(.) e 
L°°([0,r/e„],R'"), ||w„(.)||oo < 1- Then the sequence (r {x„{t / En))) ^^^^ is com- 
pact for the topology of uniform convergence on [0, T] and the limit of any converging 
sub-sequence is a solution x*{.) of the average differential inclusion (|54p . 

Proof. This proof is done as if M was M" , d the Euclidean distance and all vector 
fields had a common compact support, hence all maps a global Lipschitz constant and 
a global bound; it is left to the reader to check, as in the proof of Theorem 13.71 that 
it translates on M with any distance d described in ^2.31 

Since (Se^^g) in ((52|) can be re- written as 

^ / eg{eo + ^,x)u x 0(00 + ^,x)u 

dA V ^(eo + ^,x)+egieo + ^,x)u) ij{eo + ^,x) ' 

the norm of the difference between the right-hand sides of (Sgo.e) ^nd of 



dx _ G{9q + ^,x)u 
dA " ^^(6*0 -h^.x) 



(S«o,e) - = \: . \' A , II^II<1 (58) 



is bounded by fee for some fc > depending on the bound and Lipschitz constant of 
G and w; classical theorems on smooth dependence of solutions on the equation yield 
some c' such that, if two solutions x^{.) of (Se^^g) and Xi,{.) of (Egg g) are defined on 
the interval [0, A], with the same control A i— )■ u{X) and ^^(O) = x^{Q), then 

d{x^{\),x,{X)) <c' e for all A e [0, A] . (59) 
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Let r !->■ Xo{t) be a solution of ([51)) on [0,T]. Then A i-)- Xo(A), with 

A = T(t) = [ t:J(j;o(t))dt (60) 
Jo 

and xo{t) — xq('J'{t)) for all r, is a solution on [0, J'(r)] of 

^ G =7^£(^o) . (61) 

This is, according to ([55)1 . the average system of the fast oscillating ((55|): Theorem l3.7l 
(Point 1) yields a family of controls such that the solutions Xe{-) of (Ee^^g) with 
initial condition xq{0) and control converge to 5;o(.) uniformly on [0, J'(T)] with 
a distance less that c"e. For each e, the solution a;e(.) of (Seg e) with same initial 
condition and control satisfies ((59|) with A = J'(T); then define 



t = T(A) = - / . . (62) 

and the controls 1 1— >■ Ue{t) by Ue(A) = ?Ie(T(A)); the solutions x^{.) of (|T7|) with these 
controls are given by a?£(A) = Xs{'^{X)), and one therefore has 

d{xe{1o'P{T)),xo{T)) <{c' + c'')e, re[0,T]. (63) 

Let us now study ToCP(t). From and (15^ . one has 

dr e c^(0o + 7,2?s(A))+e5(^o + 7,2?e(A))Me(A) 

if ZU(.To(r)) 



£ \t^(6'o + 7,^o(A)) 



+ p(e,6'o,r) with A = T(t) , 



where p contains a term for the difference between this fraction with or without egue 
in the denominator, of order e by a computation similar to (|57p . and a term for the 
difference between x(A) and x{X) as an argument of ui, of order e by ([55)1 and Lipschitz 
continuity of cj; hence \p\/e is bounded by a constant independent of e. Finally, 

z:J(a;o(t))dt _ f^'^^^ dA _ ^ f[ dAd6i 



oc^(0o + ^,So(3'(t))) io c^(eo + 7,^o(A)) 27r 77 ^((?o + ^, Jo(A)) ' 

0<A<3'(t) 

and, since ZI;(a;o(.)) — '^{xvi{'y {■))), and with /i(6') is defined in ^2.4[ 

'■^aJ(xo(t))dt /"^(^^ dA 1 /■/■ dAd6' 



cJ(Jo(J'(t))) 7o ^^(^o(A)) 2^ JK<xS(r) ^^(^0+ e + ^' ^o(^)) 
1 [[ Aide 



hence the difference is less than ke for some k independent of e; the above implies that 
'XoJ'(r) — - is bounded by a constant. Finally, since Xe{.) is Lipschitz continuous 
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with constant 2£sup ||G||/fct^, one has d{x^{1 o 'J'{t)),x^{^)) < c"'e for some e and 
finally implies point 1 of the theorem, with c = c' + c" + c"' in (|62p . 

For point 2, consider (0„(.), Xn{.)^ : [0,T/e„] — )■ 5*^ x K a solution of system (|T7)) 
with e = En and some control u = Un{t). One associates to these, following ([5(1)) - 
([5^ and setting A = 3?n(i) (we write IR„ because 3? in (150)1 is constructed for system 
(Sgg^e^) and thus depends on n), solutions and control A i— >■ a;n(A) and A i— >■ Un(A) of 
(Egg^e^). The solutions A i— >■ 5;„(A) of (Egg g^) with same control and initial condition 
satisfy, according to (|59| . d(a;„(A), a;„(A)) < c'e„. By Theorem 13.71 (Point 2), the 
sequence (a;„) is compact and subsequences converge to solutions X i-^ xo{X) of (|6T)) . 

hence the same subsequences of (xji) converge as well, and, with r = Q(A) = /^^ -^g^^^^ , 
the maps r i— >■ ■z;„(Q~^(r)) = x„((Q o 3?„)~-'^(r)) converge to a solution r H> a;o(r) = 
a?o(Q~^ (''■)) of the average system (|54)) . with distance less than c"e„ for some c". Using 
the same argument as in Point 1 for ToJ'(r), one gets a bound for |(Qo3?„)~^(t) — ^| 
and, for some c'", d(x„((Q o 31„)-1(t)) , x„(-)) < c"'e„. Point 2 is proved. □ 

4.4. Dimension of £(a;). In ^3.3) and in particular in Proposition 13.101 G can 
simply be replaced with G. It is however interesting to give a more intrinsic charac- 
terization of r{9, x) and thus of dim £(.t). 

Proposition 4.4. // (117)1 and (gnj) represent the same control system, then 



dim (J2 Range ^ (61, x) 

jer 



= — 1 + rank 



{foiO,x)} U {ad}^/fe (0,x) , J e N, 1 < fc < m}) . (64) 

Proof. Straightforward computation using the fact that fo — d/d9. □ 

Note that the right-hand side is r{d,x). Proposition 13.101 applies, with this defi- 
nition of r. In particular, the "full rank case" becomes: 

Proposition 4.5. // the vector fields fo and ad-|.^/fe,r l<k<m, jeN span 
the whole tangent space of x M, then E{x) has a nonempty interior for all x. 

4.5. The function H{x,p). Instead of ([5]), H has to be taken as follows, with 
ZJ defined in (l53l): 



H{x,p) ^ uj{x) — \\{p, 

271 J Q "\ uj{e,x) 



\Ae . (65) 



The characterization of £(a;) in Proposition 13.41 is unchanged. In the "full rank case", 
the results from ij3.4l on the degree of differentiability apply without a change. 

4.6. Application to the minimum time problem. As in §3.5) but for the 

Kepler system (|47[) . let xo, xi be fixed, call the minimum time such that, from 
some 0o,6i, {9i,xi) can be reached from (^Oi^^o) in system (/Cg) (obviously Tg — )■ +oo 
as e — >■ 0) and the minimum time such that xi can be reached from xq in the 
average system ()54p . The equivalent of Theorem I3.19[ with a similar proof, using 
Theorem 1131 is: 



Theorem 4.6. In the full rank case, one has lim sup^^^Q eT^ < Tq (hence sT^ 
is bounded as s^O). If, for all e > small enough, there is a minimizing solution 
(9s, Xe) '■ [0,Te] X M and they all remain in a common compact subset of M , 

then all accumulation points (as s ^ 0) of the compact family (r — )■ Xj(|-))£>o in 
C'^([0, To], M) are minimizing for the average system and lim^^Q e — Tq . 
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The Hamiltoniaii for minimum time for the average system is given by (j65p ; one 
has to perform the time scahng described in ij4.1l to have a result Hke Theorem 13.221 
and the simple "commutation between averaging and writing Hamiltonian" noted in 
Remark 13.211 Let us translate in terms of (|47p the sufficient condition for existence 
of a Hamiltonian flow given by Theorem 13.221 

Theorem 4.7. In the full rank case, assume that {p, G{6, x)) and (p, dG/d9{9, x)) 
do not vanish simultaneously outside {p — 0}, that 9 i— )■ {p,G{9,x)) vanishes at most 
once for each {x,p) £ T*M, p ^ 0, and that ra.nkQ{9, x) = m for each {9, x) E S^x M. 
Then (|44p . with H given by (|65p . has unique solutions through any initial condition. 

The discussion of Theorem 13.221 also applies to the above; let us mention that, 
once it has been checked that, for each {x,p), {p, G{9,x)) vanishes for at most one 9, 
the other conditions are guaranteed if, (j46p holds with Q replaced by G or, in terms 
of the vector fields in (HUl) . if, for all ^ — {9, x), 



rank{/o(0, fiiO, • • • . /m (0 > ad/„ /i (C) , . . . , ad/„/,„(^)} =n+l. 



(66) 



We prove in next section that the above conditions are true at least for the planar 
control 2-body problem. 

5. Application to the controlled 2-body system. In this section we study 
some properties of the planar control system and demonstrate that it satisfies the 
condition of Theorem 13.221 on the domain of non-degenerated elliptic orbits. 

5.1. Planar control 2-body system. It is classically described by some first 
integrals of the free movement — here the semi-major axis a and the eccentricity 
vector (ex,ey) — and one angle L following the dynamics; we restrict to the set of 
non-degenerated elliptic orbits rotating in the direct sense, i.e. the state space is 
S'^ X M with M = {{a,ex,ey) G R^, a > and e^'^ + ey^ < 1}. The control 
u = {ut,Un) is expressed in the tangential- normal frame and the system reads: 



d_ 

di 







( ' ) 


( 




1 







Cy 


- „3/2 











\w{ex,ey,L)j 


V 



/ 2aEiaiex,ey,L) \ /^"A 

L) \Un) 

2 a.j, (e^; ,ey,L) hy (e^; ,ey,L) 
/ 



(67) 



with w(e2;,ej,,L) 

^ai^x, Cy, L) 
^y {^x 1 ^y : -^) 

^x{ex, ey, L) 



(1 -I- cos L + Cy sin L) 

+ e^ + 2ex cos L + 2eySmL 

y/l + + 2 ex cos L + 2ey sin L 



^yl + e'^ + 2 ex cos L + 2ey sin L 



^yl + e^ + 2 Ex COS L + 2ey sin L 



[cx + cosL) . 



[cx + cos L) . 



2 Cy + {e^ — By — 1) sin L — 2 CxCy cos L 



1 + 2ex cos L + 2ey sin L 
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^yl + e'^ + 2 Ex COS L + 2ey sin L 

2 



2ej: + (e^ -Cy + l) cosL + 2ea;ej, sinL 



1 + 2 COS L + 2 Cy s'm L 



The eccentricity e is the norm of the eccentricity vector, e — ^e^. + e^. Low thrust 
translates into ||u|| < e for a smaU e. 

Remark 5.1. This is indeed a "Kepler control system" of the type (|Tf)) except that 
a; = ^sjc?!"^ is, although strictly positive, not bounded from below by a positive con- 
stant on X M. There is such a lowerbound if one replaces M by — {(a, e^;, e^) S 
M'^, a > ande^;^ + < c} with c < 1. Strictly speaking, the results of the paper 
have to be applied in Af^, c < 1. However, Theorems 14.31 or 14.71 for instance, may 
be applied in M because each statement may ultimately be restricted to a compact 
subset of M, itself included in some M'^, c < 1. 

The Hamiltonian that both defines the average system according to ^ and yields 
the Hamiltonian system governing extremals for minimum time is given by (|65p . Since 
Jjj dL/v{ex, Gy, L) — 2n, it can be expressed as 
H{a,ex,ey,Pa,Pe,,Pey) = ^/a'H{ex,ey,apa,Pe,,Pey) with 

27r 



n{ex,ey,A,X,Y)=7r / \\{AX Y)G{ex,ey, L)\\, 



1 

2^ 

2aa/w 

G{ex,ey,L)^\ 2a.x/^ h^/^ 

2 ay/vi by/w 



5.2. Hamiltonian flow. Theorem 14.71 applies to this system. Indeed: 
Proposition 5.2. Fore each {ex,ey,a) with e^^ + Cy^ < 1 and a > 0, and each 
(A, X, Y) (0, 0, 0), the vector (AX Y)G{ex, Cy, L) vanishes for at most one angle L. 
Proof. Removing denominators, the equations Aa.a + Xsl^ + Y^Ly — and Xhx + 



Y'by = can be written: 



(2 e^A + 2(1 - e^)X) cos L + {2 CyA +2(1 - e^)Y) sin L 

= -{l + e^)A-2ex{l - e^)X -2ey{l-e^)Y 
{-2exeyX + {el ~ + 1)Y) cosL+ ((e^ - e,^ - 1)X + 2exeyY) sinL 



= 2eyX -2e,jY. 



If A = 



2exA + 2{\- e^)X 2 e^A + 2(1 - e^)^ 

-2 exCyX + {el - + i)y (^2 „ ^2 _ i)^ + 2 e^CyY 



IS nonzero, 



there is clearly at most one solution L. If A = 0, there exists A 7^ such that 

2 CxA + 2(1 - e^)X = A (-2 e^CyX + {el - + l)Y) , 
2 CyA + 2(1 - e^)Y ^ A {{el - - 1)X + 2 e^CyY) , 

and there may be a solution to the system above only if 

(1 + e^)A + 2 6:^(1 - e^)X + 2 6^(1 - e^)Y = -2 \{eyX - e^Y) 
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These three equations forms a Unear system in {A,X,Y), AI{A,X,Y)'^ — with 

/ 2e^ 2{l ~ + Xe^cy) -A(e2 - + 1) \ 
M = 2ey -\(el - - 1) 2(1 - - Xe^Cy) 

V (l + e2) 2{e,{l~e') + Xey) 2 (ey(l - e^) - Ae,) / 

A brief computation gives detM = (1 — e)"^(l + e)'^(A^ + 4), strictly positive when 
< e < 1. Hence M{A, X, F)^ = impHes {A, X, F) = 0. □ 
Since the rank of G is obviously equal to 2 and the rank of {G, dG/dL} equal to 3 
for any {ex, Cy, L), the hypotheses of Theorem 14 .71 are satisfied by the planar control 
2-body system, and it guarantees existence of a flow for the Hamiltonian system 
governing the extremals of minimum time for its average system. 

6. Conclusion. We have defined an average control system in a way that differs 
from the usual use of averaging in control theory, and given its basic properties: an 
equivalent for control systems of the averaging theorems for ODEs that justifies the 
average control system as such, and statements on its description and regularity. It 
has natural applications to approximation of minimum time control. It has already 
allowed us to give (with restrictions on the eccentricities coming from the point raised 
in Remark 15. ip a proof that the minimum time between 2 ellipses grows like 1 /e for 
the planar 2-body problem [B]. 

The present results are however mostly a starting point. The regularity of H has 
to be further explored when the conditions of Theorem 13. 221 do not hold, see the last 
paragraph of ^ Explicit computation of the average system and its extremals for 
the 2-body problem has to be conducted. 
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and the two anonymous referees from this journal for extremely constructive reviews 
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Appendix. Proof of Proposition 13.161 We first state two lemmas. 

Lemma A.l. Assume that X Cz Z and (j39p is satisfied. There is a neighborhood 
U of X in O'^ and a smooth map x '■ U such that, for {0,X) G U, one has 

\l{e,X) = only if 9 = x{X), and 



^mnx) 



ymx),x) 



0, X eU. 



(A.l) 



Proof From ([Sla), Z = {(6'^X) G x O'^ , V(6',X) = 0} is a smooth submani- 
fold of X O'^ and from (|39lc). Z given by (|34| a smooth submanifold of C^, both of 
dimension d+1 — m, and the projection vr : S*^ x O'' — > O'^ induces a diffeomorphism 
Z ^ Z whose inverse is of the form x i— > (x(a;),a;) with x ^ smooth map Z — > 5^ 
that satisfies, for all X € Z: V(6', A) = if an only if 6* = x{x)- 

Consider the map T : x O'^ ~> R given by T{e,X) = {%{0,X) |V(6I, A)). Let 
X be in_Z; _since V(x(A), A) = 0, one has T(x(A), A) = and dT/de{x{X),X) = 
II ^(x(A), A)|p, nonzero from assumption (|39lb): the implicit function theorem 
yields a unique map x from a neighborhood U of A in O'^ to a neighborhood of 
x{X) in such that 9 — £(A) solves T{9,X) = 0; it must therefore coincide with x 
in U O Z and satisfies the lemma. □ 



AVERAGED CONTROL SYSTEM 



23 



Lemma A. 2. Assume that X £ Z and (|39p is satisfied. There exist a neighbor- 
hood U of X in O'^ , local coordinates Xi, . . . , defined on U , and smooth maps 

( xi \ 

P: U ^ SO{m), a: U ^R, and W: X U such that, with Xi= : 

\Xrn-l/ 

(A.2) 
(A.3) 
(A.4) 



yi9,X)^P{X) 
^P{X) 



Xi 

a{X) {9-x{X)) 

Xi 

^ 

Om-1 



-{e-x{x)r w{0,x) 

iO-x{X)) Wi{0,X) 



with w^{e,x) = ( ^'j-^^ ) + {e-x{x)) w{e,x) 



in X U, where a is bounded from below: < ao < o:{X), X ^U. Furthermore, for 
a constant > 0, one has, for all {6, X) E x U , 



and 



M0,x)\\ > i^sV ll^ilP + " W - x{x)r , 

Xi^O ^ \\Wi{9,X)\\ >K3. 



(A.5) 
(A.6) 



Proof. The map X t-s- ^{x{X), X) is nonzero for X — X, hence it does not 
vanish on a sufHciently small neighborhood U of X , and one may write 

^mX),X) = PiX) ) , a{X) > ao > . (A.7) 

Define vi, . . . ,Vm, smooth maps S*^ x [/ M by 

/ v^{e,x) \ 

: \=p-\X)y{9,X). (A.8) 

V v^{e,x) J 

For i between 1 and m - f, ^{x{X),X) ^ from ((XT)) . and w,(x(X),X) = 
from Lemma I A. II and. using ipSl a). the rank of the map X i-> (wi(x(^), ^), • . • , 
Vm~\{xXX^,X')') is m — 1 at X = X: on a possibly smaller neighborhood U , there 
are local coordinates x\, . . . ,Xd such that Vi{6,X) = Xi + {6 ~ x{X)) Wi{d,X) for 
i < m — 1 and for some smooth Wi; substituting (IA.7I) and (jA.SP in (lA.ip implies 
v^{x{X),X) =. 0, hence «„(0, X) = a{X) [9 ~ x{X)) + W^{e, X) {B - x{X)f for a 
smooth M/^m; (|A.2p is proved. 

Possibly restricting [/ to a subset with compact closure, || VF((?, A")|| is bounded 
on 5^ X [/; if 16* - x{X)\ < i^o/ max ||W^||, then (jA.sp holds with K3 = ^ according 
to (jA.2p ; on the set where \9 — x{X)\ > ^ao/max||H^||, V does not vanish and hence 
(||Xi||2 + a(X)2 (6* - x{X)f)^/^/\\y{0,X)\\ is bounded from below; ^EM is proved, 
with K3 smaller than this bound and than i. From (jA.4p . Wi{x{X), X) 7^ because 
a does not vanish; from assumption (|39lb) and (jA.3P (where Xj — if X — X), 
Wi{d,X) ^ if 7^ xiX), hence Wi does not vanish on x {X}; it is therefore 
bounded from below on x U with U a small enough neighborhood of X: (|A.6|) 
holds with K3 smaller than this bound. □ 
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Proof of Provosition \3.l6[ We use [— 7r,7r] instead of [0,27r] as an interval of 
integration. Let h E S.'^ , with ||/i|| = 1. From one has, for some constant K 

using bounds on the derivatives of the smooth V, 



|dH(X)./i- dH(y)./i| < 



1 

2^ 



/5V 



lx).h^--{e,Y).h 



1 

2^ 



lY).h 



9V 

ax' 
v(^^,x) 



V(0,X) 



M9,X)\ 
y{0,Y) 

\\y{e,x)\\ \M9,Y)\\ 



<K\IX-Y\\ + ^ 



" y{0,x) 

-J|V(0,^)l 



d6i - 



V(6i,F) 

-J|v(0,r)|' 



de 
de* 



FinaUy, defining 

\/{cp,X)=y{x{X)+ip,X), Wi{ip,X)^W,ixiX) + ^,X) , (A.9) 
and making a different change of variables in the last two integrals, one has 

K r,, y{ip,x) v(^,y) 



\dH{X).h - dH{Y).h\\ < K\\X - Y\\ 



27r 

<K\\X-Y\\ + - 

TT 



'||v(^,x)|| |iv(^,y)| 
||v(^,x)-v(^,y)| 



di^ 



|V(^,X)|| 



d^ (A.IO) 



where the last inequality uses the fact ||||^ — ||^|| < 2 min{ , }, and also 

holds with ||V((/?, instead of ||V((/3, in the denominator. Now let us use 
Lemma IA.21 let X — {xi, . . . , Xd) and Y = (yi, . . . , yj) in these coordinates; from 
()A.3D . one has, with Wi defined by ()A.9|) . 

yi^,X) ^ PiX)[(^A + ipWi{ip,X)], yiif,Y) = PiY)[(^A + ^Wiiip^Y)]. (A.ll) 



Hence V(^, X) - V(^, Y) = {P{X) - P{Y))P{X)-^y{^, X) 

+P{Y)[^{Wi{^,X)-Wi{^,Y)^ + 



and finally 

||V(^,X)-V(y.,y)H 
l|V(^,^)ll 



X^-Yi 




<\\P{X)-P{Y)\\ 



\^\\\W^{^,X)-W^{ip,Y)\\ 

\fj{v.x)\\ 



\Xi-Yi\ 



\\'J{v.x)\\ 

(A.12) 

Two cases are to be distinguished: 

(i) If Xi = Yi = 0, then </? factors out of \/{ip, X) and Y) in (jA.lip and the 
last term in (|A.12p is zero: according to (jA.6l) . the integrand in (|A.10p is bounded by 



\\P{X)-P{Y)\\ + 



\Wi{ip,X)-W^{ip,Y)\\ 

K3 



and finally \dH{X).h - dH{Y).h\ < K \\X - Y\\ with a constant K that depends only 
on V, the open set U and the coordinates. 

(ii) If Xi 7^ (or Yi ^ 0, interchanging X and F), then (|Al2| . using (jA.5|) . 
implies that the integrand in (jA.10| is bounded by 
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but the same is also true replacing a{X) with a{Y) and ||Xi|| with UliU ; hence, since 

||a-6|p < 4max{||a||2, the last term maybe replaced by ^ || Xiiy|f^+£oy^ ' 

whose integral between — tt and tt is equal to 



which is less than — liH (fci + A;2 In pj^^ry^) for some fci, ^2 when, say, ^^^'^J^'^^ < 1, 
and finally, since \\Xi — Yi\\ is less than \\X — Y\\ and u uln(l/w) is nondecreasing, 
less than ||X - F||(fci + k2 In ||x-y|| )- 

Cases (i) and (ii) do imply (|40ll . possibly restricting U so that In ux-y\\ '-' 
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